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ABSTRACT 

Cool-core clusters (e.g., Perseus or M87) often possess a network of bright gaseous fila- 
ments, observed in radio, IR, optical and X-ray bands. We propose that these filaments 
are powered by the reconnection of the magnetic field in the wakes of buoyant bubbles. 
AGN-inflated bubbles of relativistic plasma rise buoyantly in the cluster atmosphere, 
stretching and amplifying the field in the wake to values of /? = SirPgas/B'^ ^ 1. The 
field lines in the wake have opposite directions and are forced together as the bubble 
motion stretches the filament. This setup bears strong similarity to the coronal loops 
on the Sun or the Earth magneto-tail. The reconnection process naturally explains 
both the required level of local dissipation rate in filaments and the overall luminosity 
of filaments. The original source of power for the filaments is the potential energy of 
buoyant bubbles, inflated by the central AGN. 

Key words: methods: numerical - galaxies: clusters: intracluster medium - X-rays: 
galaxies: clusters 



1 INTRODUCTION 

Networks of bright gaseous filaments are ubiquitous in the 
centers of cool-core clusters (e.g., McDonald et al.|2010 |. Hq 
filaments around NGC1275 in the Perseus cluster are per- 
haps the most famous example (e.g., |Minkowski|1957|[Lynds 
|1970 |. These filaments are observed in many bands/lines, in- 
cluding C O (e.g.,|Lazareff et al.]|1989| [Salome et a"L]|2006[), 



NIR lines ([Mittalet al.|2012[ ), optical lines (e .g.,|Conselice^ 



Gallagher, & Wyse 20011 and soft X-rays (Fabian et al 



2003), suggesting a multi-temperature gas sharing approxi- 



mately the same space within the cluster. For a recent sum- 
mary of observational results on NGC1275 and M87 fila- 
ments see 



e.g.. 



Fabian et al. (20111; Werner et al. (20121 



and references therein. Below we discuss NGC1275 and M87 
collectively, under implicit assumption that the same mech- 
anism is responsible for the filamentary structures in both 
objects (and also in other cool-core clusters). 

The source of energy powering the filaments is a long- 
standing problem. The bolometric luminosity of the fila- 
ments in the NIR-optical band could be at the level of 10- 
20% of the total X-ray luminosity of the cluster core. Vari- 



Bregman, fc Seabl|1988l, photoionization by optical/UV or 



X-ray radiation (e.g.,|Heckman et al.|1989[|Voit, Donahue, fc" 



Slavin 19941, thermal conduction (e.g., Boehringer & Fabian 



19891. Ferland et al. (20091 argued that the spectra of the 



ous scenarios have been considered, including shocks ( David, 



outer filaments require the line excitation by energetic par- 
ticles, although not all line ratios are consistent with this 
scen ario ([Mittal et ar||2012[ ) . Recently [Fabian et al] ( |2011 1 
and [Werner et al.| ( |2012[ ) suggested that the filaments are 
powered by the hot ICM, which penetrates into the fila- 
ments via turbulent reconnection. 

The filaments are long and thin, probably consisting of 
many threads ( [Forman et al.|2007[ [Fabian et al.|2008[ ). This 
suggests that the magnetic field is playing a role. The role 
of magnetic fields and in particular reconnection as a source 
of energy for filaments has been considered in, e.g., |Soker fc| 
[Sarazin[ ( [1990[ ) and [Jafelice fc Friaca[ ( [1996| ). It was assumed 
that an inflow of cooling gas (in the frame of original cooling 
flow model) increases the magnetic energy density in the core 
of the cluster. The relative contribution of the magnetic field 
to the energy density is further amplified by the radiative 
cooling losses of the gas thermal energy. 

Here we consider a different scenario, in which buoyant 
bubbles of relativistic plasma stretch the magnetic field lines 
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and drive the fields of opposite direction together. In this 
model the AGN-inflated bubbles provide the energy that 
powers the filaments. A schematic picture of this process is 
shown in Fig. [T] 

The structure of this paper is as follows. In !|2]we briefiy 
summarize the relevant properties of AGN-inflated bubbles. 
In !|3] we discuss the amplification of the magnetic field by 
the rising bubbles. In iplwe provide an order-of-magnitude 
estimate of the rate of energy dissipation by reconnecting 
magnetic fields in the bubble's tail and the resulting lumi- 
nosity of the filaments. In ^we discuss the overall energetics 
of filaments and the most basic properties of our model. Our 
findings are summarized in !|6] 



2 BUOYANT BUBBLES 

Observations suggest that AGN activity regulates the ther- 
mal state of the gas by injecting energy into the intra-cluster 
medium in the cores of relaxed clusters, where radiative cool- 
ing time is often as short as few 10* years. Bubbles of rel- 
ativistic plasma are inflated by a supermassive black hole 
and rise buoyantly through the gaseous atmosphere, leading 
to a number of spectacular phenomena such as expanding 
shocks. X-ray-dim and radio-bright cavities, X-ray-dim and 
radio-dim "ghost" cavities (aged versions of "normal" cav- 
ities) and filaments of cool gas in the wakes of the rising 
bubbles formed by the entrained low-entropy material from 
the core ( [Churazov et al.||2000l |2001| . With Chandra and 
XMM-Newton, these features are now studied in great de- 
tail in many systems. 

Observations further suggest that large fraction of en- 
ergy output of the AGN goes into the enthalpy of the bubbles 
H = 7i,PVi,/(76 — 1), rather than into shocks. Here 75 is the 
adiabatic index of the gas inside the bubble (75 = 4/3 or 5/3 
depending on whether a relativistic or a non-relativistic gas 
is considered), P is the pressure inside the bubble and Vb 
is the bubble volume. The partitioning of the AGN energy 
between shocks and bubble enthalpy depends on the energy 
injection rate, duration of the AGN outburst and initial con- 
ditions (e.g., Forman et al.|2007 20131, but fiducial models 
predict that about 70% (and certainly more than 50%) goes 
into the enthalpy of bubbles. Furthermore, the lack of very 
strong shocks around observed bubbles suggests that the 
thermal gas pressure of the ICM supporting the bubble, can 
be used in the above expression for the enthalpy H (i.e., 
P ~ Pgas, the bubble and the ICM are in pressure balance). 

The bubbles then serve as a reservoir of potential of 
energy ~ H, deposited by the AGN. The dynamics of the 
bubble rise is set by the competition of the buoyancy force 
and the drag from the ambient gas. Even if we consider only 
the hydrodynamic drag (i.e., ignoring a possible contribu- 
tion of magnetic fields) the rise velocity is expected to be 
subsonic |Churazov et al.|2000[ ). Indeed, the buoyancy force 
is Fb ~ pVhQ, where g is the gravitational acceleration, and 
it is balanced by the ram pressure of the ICM (inertial drag 
force) Fram ~ Apvt^^^, where vterm is the bubble's termi- 
nal velocity, A is the cross section of the bubble. Equating 
Fram and Fb gives an expression of the terminal velocity 
Vterm ~ V 9 R, where R is the bubble radius. The terminal 



velocity should be subsonic/transsonic as long as the bub- 
ble radius does not exceed the pressure scale height of the 
atmosphere. Assuming that the bubble is moving subsoni- 
cally and does not mix with the ambient ICM, the volume 

f P\~^ 
of the bubble expands adiabatically Vt = Vtfi 



where P = Pir), Po = P{ro) is the ICM pressure and ro 
and r are the initial and current distance of the bubble from 
the cluster center. For simplicity, we assume below a power- 



law dependence of pressure on the radius P = Po 



ro 



Typically a ~ 0.7 — 1 for the relevant range of radii in cool 
core clusters. For example, using radial density and temper- 



ature profiles from Churazov et al. ( 2003 1 and ( Forman et 
al.|2007l we obtained a = 0.8 and 0.9 for the Perseus cluster 



and M87 respectively. 

The ambient material and the bubble itself can be 
threaded by the magnetic fields. As the bubble rises, the 
magnetic field is amplified - this process is discussed in the 
next section. 



3 RISE OF THE BUBBLE 

The role of magnetic fields in the evolution of buoyant bub- 



bles has been considered in, e.g., Ruszkowski et al. (2007 



2008| ; lO'NeiU, De Young, fc Jones] ( |2009| ). Here we concen 
trate specifically on the threads of the magnetic field in the 
wake of a rising bubble. 

A sketch of the configuration is shown in Figjl] We as- 
sume that the bubble advects a lump of the ICM threaded 
by magnetic field lines, which are anchored to the gas in the 
cluster core. We assume that initially the reconnection of the 
magnetic field can be neglected. As the bubbles rise, the ad- 
vected fiuid elements and the magnetic field frozen into them 
are stretched by the bubble motion. We start by considering 
the evolution of such an advected fiuid element occupying 
a volume 1/ < 14. As the fluid element moves from ro to 

r, its volume expands adiabatically V = Vo [ jr] '' , where 
7 is the adiabatic index of the ICM. The linear size of the 
stretched fluid element along the direction of motion can be 
estimated as Z ~ 7?o -I- r — ro, where i?o is the initial size 
of the fluid element. In the limit r ^ ro, Z ~ r. The cross 
section A of the fluid element in the perpendicular direction 
is then 



A' 



V_ 
1 



I 



Vo 



ro 



^-1 
oc ri 



(1) 



Thus, for a < 7, the cross section of the fiuid element shrinks 
as the bubble rises. 

The stretching of the fiuid element will align and am- 
plify the magnetic field B. From the conservation of the 
magnetic flux through the cross section of the advected fluid 
element: 



B_ 

B'o 



A 



J_ 
Ro 



r 
R'o 



ro 



(2) 



where Bo is the initial magnetic field and Ao ~ Vo/ Ro the 
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Figure 1. Qualitative picture of the current sheets formation by AGN-infiated buoyant bubbles of relativistic plasma, rising in the cluster 
atmosphere. Nearby elliptical galaxy M87/Virgo is used in this example. Left: Morphology of soft X-ray filaments in M87 ( Form an et al.| 
|2007| | and overall morphology of the radio emitting plasma jOwen, Eilek, fc K assini 2000), superposed as contours. Optical filaments are 
largely co-spatial with X-ray filaments. Buoyant bubbles rise in the atmosphere, entraining the low entropy gas from the core (|Chura20v| 
|et al .'2000 2001 ) and stretching/squeezing the fluid elements in the wake. Radio emission traces the distribution of the relativistic plasma 
produced by the AGN. Right: Schematic evolution of the magnetic field in the wake. As the bubbles rise they stretch the magnetic field 
lines in the entrained fluid elements, thus increasing the strength of the field. The field lines, anchored to the gas in the cluster core, have 
opposite directions in the wake. They are forced together as the bubble rises. This setup bears strong similarity to the coronal loops on 
the Sun or the Earth magneto-tail, where reconnection is taking place. 



initial cross-section. The corresponding magnetic energy is 



B 



Bl 



E^^^V=^Vo 



P 



Stt 



where /3o 



TT "Rl \Po 



PoVo r^ 
Po Rl \ro 



(3) 



Stt-^ ~ 100 (e.g., 

^0 



Carilli & Taylor 



20021 is 

the /3 parameter of the ICM near the initial position ot the 
bubble. Using this expression, we can estimate the maximum 
distance rmax from the cluster center the bubble can reach 
- this is the radius where the buoyancy force Ft ~ pVtg ~ 

Vb —r- ~ Vb — is equal to Fl ~ — ; — ~ . In the limit of 

dr r dr r 

rmax 3> ro,Ro, the equality Ft ~ F^ is reached at 



ro /3o 



14,0 Ro 



2 + 



:v/7-a/7^ 



(4) 



Vo rl 

At this radius, the value of /3 in the stretched fluid element 
is 

/3(.™.)~^~/3o^(^/?o^^J . (5) 

Setting Vtfi ~ ''o (i-e., initial bubble size is comparable with 
the initial distance from the cluster center) and Vb ~ i?0) 
we get 



P [rmax) 



ro (A, 
/3o (Po 



rp 
Ro 

rp 

Rp 



^/-i-o'f-ih 



2-, + a(-,-l--,/76) Jj2 



For a set of fiducial values a 
A/2,, Po ~ 100, we get 

rmax = 10 ro (^ j 

Pirmax) - (j^y . 
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0.85, 7 = 5/3, 76 



(6) 
(7) 

(8) 
(9) 



Thus, it is reasonable to expect the bubble to rise a distance 
of order r ~ 10 ro before the buoyancy and magnetic forces 
come into balance. At this point, the /3 parameter in the 
stretched fluid elements approaches unity. Circumstantial 
evidence for the magnetic field energy density comparable 
to the ICM thermal pressure was indeed presented (based 
on a different argument) in Fabian et al. (2011 1; Werner et 



4 RECONNECTION IN THE FILAMENTS 

Once the bubble is at rmax, further stretching of the field 
lines is not possible. The bubble (or fiuid elements attached 
to it) would "hang" on the magnetic field lines. However, 
the field lines in the stretched fluid elements will have op- 
posite directions and are forced together by the shrinking 
cross section of the fllament. The conflguration bears strong 
similarity to the solar coronal loops (e.g 



Kopp fc Pneuman 
19761 or the Earth magneto-tail (e.g., Nishida||2000l, mak- 



ing the fllament prone to reconnection. As the anti-parallel 
fleld lines come together, current sheets are formed, with 
an inflow of magnetic energy, which is eventually dissipated 
there. The release of magnetic energy allows the bubble to 
rise further. 

Magnetic reconnection in both coUisional (MHD) and 
collisionless plasmas proceeds at a rate that is independent 



of Ohmic resistivity (Uzdensky, Loureiro, & Schekochihin 



20101 or other aspects of plasma microphysics (Rogers et 
al. 2001). Namely, one can write a rough estimate of the 
magnetic energy inflow per unit surface of the reconnecting 
layer as follows: 



; evA 



?1 
87r' 



(10) 
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where va = \l is the Alfven speed, n is the gas 

47rn/imp 

particle density, /i mean particle atomic weight and e is a 

dimensionless coefficient varying between e ~ 0.01 for col- 

lisional plasmas (Uzdensky, Loureiro, fc Schekochihin|2010 



Loureiro et al.|2012| Bhattacharjee et al.|2009[|Daughton et| 



al. 2009[ Loureiro et al. 20091 and e ~ 0.1 for collisionless 
ones ( Birn et al.|2001 1. Assuming that /3 ~ 1, we can replace 
the magnetic energy density with the thermal energy density 

B^ 

— ~ nkT and the Alfven velocity va with the sound speed 

8tt 



Ij—^ — . This gives an order-of- magnitude estimate of 
the surface influx of energjr] 



: eCsUkT. 



(11) 



For the NGC1275 and M87 the estimates of the total emitted 



surface flux by the filaments are available ( [Fabian et al.|20lT 
Werner etal. |2012[ ): Le, " ' 



10"^ ergs s"' cm"^ ~ 0.2csnkT 
^ ~ O.lCsTifcT, respectively. 



and ~ 2.2 10"'' ergs s"' cm 
Thus, there is an interesting order-of-magnitude agreement 
between the amount of energy that can be produced by fast 
reconnection and the amount of energy emitted by the fil- 
aments, i.e. — ^^ « — — . If the reconnection rate is on the 

^em 0.1 

stronger side of the possible values, viz. e ~ 0.1, the energy 
release from reconnection is comparable to the cooling losses 
of the filaments. 



5 DISCUSSION 

The overall energetics of the cool cores are believed to be 
determined by the balance of the AGN activity and gas 
cooling. In other words, one can assume that the cooling 
losses are approximately matched by the amount of mechan- 
ical energy pumped by the AGN into the gas in the form 
of relativistic bubbles. Observations suggest that a signifi- 
cant (if not dominant) fraction of the AGN energy goes into 
the enthalpy of the bubbles rather than in the shocks (e.g., 
Churazov et al.||2002 L This means that potential energy of 
underdense bubbles created by the AGN per unit time ap- 
proximately matches gas cooling losses. The estimates in ^ 
suggest that by the time when (3 reaches 1, the bubble has 
moved to r^ax ~ lOro. Let us estimate the ratio /s of the 
magnetic energy of the stretched fiuid element at this mo- 



^ Note that while the dissipation rate of the magnetic energy in 
the reconnection process does not have to be the same as the 
reconnection rate, it is reasonable to expect that they are com- 
parable (there is some numerical evidence in support of this, e.g., 
[Loureiro et al.^^2012j . In our simple estimates we have absorbed 
both the reconnection and the dissipation rate into the e param- 
eter. 



ment to the initial enthalpy of the bubble Ho = -^^PoVt^o'- 



fB = 



Eb 



lb 



PV 



lb P{rmax) PoVbfi 



lb 



0.1, 



(12) 



7b — t i / r,nax 

lb P{rmax) V '^0 

using eq.( [3|7|8"| ) and neglecting dependence on Ro/ro. At 
rmax, the fraction of the remaining enthalpy of the bubble 
is 



h 



H 



PVb 



Ho PoVi,,o 



ro 



^b 



0.5. 



(13) 



The rest of the initial enthalpy has already been transferred 
to the gas via hydrodynamic drag, potential energy of the 
uplifted gas, magnetic energy, excitation of y-modes, which 
then dissipate in the ICM ( [Churazov et al.[[2002[ ). Compar- 
ison of /h and /b suggests that by the time the bubble 
reaches Vmax, about 20% of its available energy will have 
gone into magnetic energy forced into its tail. The luminos- 
ity of the filaments from NIR to optical bands amounts to 
10-20% of the bolometric luminosity of the cluster cores. 
This means that 10-20% of the potential energy available 
conversion into magnetic energy should indeed go into re- 
connection and the associated heating. When the reconnec- 
tion releases magnetic energy, the bubble continues to rise 
beyond rmax- The partitioning of the remaining energy de- 
pends on the configuration, but it is likely to stay at the 
same level. 

Note that it very likely that in real clusters there is 
a considerable spread in the values of initial parameters, 
such as, e.g., /3o, ro/Ro- This suggests a large spread in the 
appearance of the filaments in different clusters or even of 
filaments/bubbles in the same cluster. 

In our simple scenario, only the regions where the con- 
figuration of magnetic field is favorable for reconnection are 
observed as filaments. The energy of the magnetic field is 
the principal source of energy for the outer filaments, rather 
than ICM thermal energy or photoionization, although both 
can contribute. The implication is that the filaments do not 
necessarily grow in mass with time, instead they thermalize 
and emit the bubble energy, mediated by magnetic fields. 

The cooling of the gas is not the central element of our 
model (cf. Soker & Sarazin 1990 Jafelice fc Friaca 1996]), 



in the sense that the main driver of the reconnection is the 
stretching of the field lines by the bubbles, rather then the 
loss of thermal energy by the cooling gas. It is neverthe- 
less clear that a large amount of cool gas is present in these 
systems. In the simplest scenario, this gas intercepts, ther- 
malizes and re-emits the released energy of the magnetic 
field. The discussion of how the magnetic energy is split be- 
tween thee kinetic energy of the gas, its thermal energy and 
non-thermal particles, and of actual excitation of the opti- 
cal lines is beyond the scope of this letter. We nevertheless 
note that the presence of non-thermal particles may help 



explain many properties of the emission spectra ( Ferland et 
|aH2009| . 

As a speculative extension of this qualitative model, we 
note that the gas leaving the reconnection region will have 
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velocities of order va, which is also of the order of the sound 
speed. The outflow is typically bi-directional, i.e., along the 
filaments. Reconnection in extended current sheets is typ- 
ically accompanied by generation and ejection of copious 
number of plasmoids (Loureiro et al.|2012 Bhattacharjee et 



al.|2009|[Daughton et al.|2009[|Samtaney et al. '2009), some 



of them very large (Loureiro et al. 20121. If the filament 
is aligned along the line of sight towards an observer, this 
may lead to the appearance of gas lumps moving towards 
the core, away from the observer with a speed that can be 
as large as ~ 10^ km s~^ for the hot gas. There is a so- 
called High Velocity system (HV) in the core of the Perseus 
cluster - a line emitting region in the core of NGC1275, 
with the recession velocity ~3000 km/s larger than the sys- 
temic velocity of NGC1275 (e.g., )Min kowski 19 57|), which is 



is nevertheless located in front of NGC1275 (e.g., De Young, 



Roberts, & Saslaw 1973 1 and therefore is moving towards 



the nucleus. While the infalling velocity of HV seems to be 
too large for a conceivable plasmoid-ejection mechanism, it 
is nevertheless interesting to note that in some favourable 
configurations, high-velocity gas lumps can be observed. 



6 CONCLUSIONS 

We argue that buoyant bubbles in the cores of galaxy clus- 
ters stretch the fluid elements advected from the core, form- 
ing gaseous fllaments and aligning and amplifying the mag- 
netic fleld in these filaments. The field grows to /? ~ 1 after 
the bubbles rise a distance of order 10 times their size. The 
fleld lines in the wake of the bubble are anti-parallel and are 
forced together. This setup bears strong similarity to the 
coronal loops on the Sun or the Earth magneto-tail. The re- 
connection process can naturally explain both the required 
level of local dissipation rate in fllaments and overall energy 
balance. In this model, the ultimate source of power for the 
fllaments is the potential energy of buoyant bubbles, inflated 
by the central AGN. Of order 10% of the total mechanical 
energy deposited by the AGN in the form of relativistic bub- 
bles can be converted into the emission from the filaments. 
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